AD-A227  205 


OWC  F/LE  COPY 

NASA  Contractor  Report  182003 
ICASE  Report  No.  90-19 


(D 


ICASE 


SURFACE  —  COOLING  EFFECTS  ON  COMPRESSIBLE 
BOUNDARY  —  LAYER  INSTABILITY 


S.  O.  Seddougui 
R.  I.  Bowles 
F.  T.  Smith 


DTIC 

ELECTE I 
0CTO3 1990 


Contract  No.  NAS 1-1 8605 
February  1990 


Institute  for  Computer  Applications  in  Science  and  Engineenrtg 
NASA  Langley  Research  Center 
Hampton,  Virginia  23665-5225 

Operated  by  the  Universities  Space  Research  Association 


fUASA 

National  Aeronautics  and 
Space  Administration 

Langley  Research  Center 

Hampton,  Virginia  23665-5225 


"distribution  statement~a" 

Approved  for  public  release; 

D:°.trib,i'iioa  Unlimited 


90  -» 


r 


Surface  -  Cooling  Effects  on  Compressible 
Boundary  -  Layer  Instability 


S.  O.  Seddougui1 


Institute  for  Computer  Applications  in  Science  and  Engineering  ^7 
NASA  Langley  Research  Center 
Hampton,  VA  23665 


R.  I.  Bowles  and  F.  T.  Smith 
University  College  London 
Gower  Street 
London  WC1E  6BT 


Abstract 


S  Accession 

For 

:  i 

'1 

•1  r- 

L 

i  Oil _ 

:  r.  '  ' 

/ 

' t 

l.y  Code 

"•  /  r-  >•  - 

.....  w  j 

;  1 

;  mi/or 
•: :T 

A'L 

i 

..i _ 

The  influence  of  surface  cooling  on  compressible  boundary-layer  instability  is  dis¬ 
cussed  theoretically  for  both  viscous  and  inviscid  modes,  at  high  Reynolds  numbers. 
The  cooling  enhances  the  surface  heat  transfer  and  shear  stress,  creating  a  high- 
heat-transfer  sublayer.  This  has  the  effect  of  distorting  and  accentuating  the  viscous 
Tollmien-Schlichting  modes  to  such  an  extent  that  their  spatial  growth  rates  become 
comparable  with,  and  can  even  exceed ,  the  growth  rates  of  inviscid  modes,  including 
those  found  previously.  This  is  for  moderate  cooling,  and  it  applies  at  any  Mach  num¬ 
ber.  In  addition,  the  moderate  cooling  destabilizes  otherwise  stable  viscous  or  inviscid 
modes,  in  particular  triggering  outward-traveling  waves  at  the  edge  of  the  boundary 
layer  in  the  supersonic  regime.  Severe  cooling  is  also  discussed  as  it  brings  compressible 
dynamics  directly  into  play  within  the  viscous  sublayer.  All  the  new  cooled  modes  found 
involve  the  heat- transfer  sublayer  quite  actively,  and  they  are  often  multi-structured  in 
form  and  may  be  distinct  from  those  observed  in  previous  computational  and  experi¬ 
mental  investigations.  The  corresponding  nonlinear  processes  are  also  pointed  out  with 
regard  to  transition  in  the  cooled  compressible  boundary  layer.  Finally,  comparisons 
with  Lysenko  and  Maslov’s  (1984)  experiments  on  surface  cooling  are  presented. 


’This  research  was  supported  by  the  National  Aeronautics  and  Space  Administration  under  NASA  Con¬ 
tract  No.  NAS1-18605  while  the  authors  were  in  residence  at  the  Institute  for  Computer  Applications  in 
Science  and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23665. 
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1.  INTRODUCTION 


Surface  cooling  is  fairly  often  applied  on  supersonic-hypersonic  flight  vehicles,  and  in 
some  cases  it  is  essential  for  the  preservation  of  the  vehicle  itself,  especially  at  higher  Mach 
numbers.  The  effects  of  such  cooling  on  the  instability  and  transition  properties  of  the  local 
compressible  boundary-layer  flow  are  then  of  much  concern.  Given  the  practical  importance, 
it  is  perhaps  surprising  that  there  are  not  more  experimental,  theoretical  and  computational 
studies  of  surface-cooling  effects  on  instability  and  transition;  see  the  experiments  of  Van 
Driest  and  McCauley  (1957),  Wisniewski  and  Jack  (1960),  Richards  and  Stollery  (1966), 
Stetson  and  Rushton  (1967),  Van  Driest  and  Blumer  (1968),  Kendall  (1975),  Demetriades 
(1978),  Lysenko  and  Maslov  (1981,1984),  and  the  theoretical  suggestions  and  computations 
of  Lees  and  Lin  (see  Lin  1955),  Reshotko  (1963),  Mack  (1969,1975,1984,1986),  Maslov  (1974), 
Gapanov  and  Maslov  (1980),  Malik  (1987).  The  flow  phenomena  produced  can  be  complex, 
and  the  findings  and  suggestions  from  the  works  above  are  a  little  mixed,  but  Lysenko  and 
Maslov  (1984)  seem  to  provide  a  reasonable  overall  summary,  namely  that  the  so-called  first 
modes  tend  to  be  stabilized  by  surface  cooling  whereas  the  second  modes  are  destabilized  to 
a  limited  extent.  Also,  islands  of  first-mode  instability  have  been  predicted  by  parallel-flow 
theory  at  the  lower  Reynolds  numbers,  although  these  suggestions  are  questionable  due  to 
the  nonparallel- flow  effects  present  then:  see  also  later. 

The  present  theoretical  investigation,  in  comparison,  suggests  that  the  effects  of  surface 
cooling  can  be  much  more  pronounced  and  “dangerous.”  In  particular  it  is  found  that  the 
cooling,  even  if  only  moderate,  so  distorts  the  viscous-inviscid  Tollmien-Schlichting  (TS) 
modes  [first  modes]  that  their  spatial  growth  rates  can  become  comparable  with  or  larger 
than  those  of  the  inviscid  modes  [second  and  higher  modes]  described  above,  at  high  Reynolds 
numbers.  This  reversal  of  the  usual  ordering  of  the  instabilities  arises  because  of  the  enhanced 
heat  transfers  and  shear  stresses  produced  in  the  basic  boundary-layer  flow  at  the  surface 
by  the  presence  of  cooling,  which  results  in  a  sublayer  of  high  heat  transfer  emerging  on  the 
surface.  The  increased  surface  shear  stresses  cause  the  typical  viscous-inviscid  TS  wavelength 
to  be  decreased  and  the  flow  structure  associated  with  the  TS  wave  is  soon  altered,  as  the 
surface  is  cooled,  bringing  about  the  rise  in  growth  rates,  as  well  as  rendering  previously 
stable  modes  unstable. 

The  present  predictions  appear  to  agree  broadly  with  Lysenko  and  Maslov’s  (1984)  sum¬ 
mary  finding  above,  and  with  their  experiments  for  supersonic  boundary  layers  (see  also 
the  comparisons  in  §6  below),  in  that  the  wavenumbers  and  frequencies  associated  with  the 
boundary-layer’s  first-mode  instability  increase  with  surface  cooling,  at  a  given  Reynolds 
number.  There  is  however  some  possible  disagreement  as  well  (which  may  be  a  matter  of 
interpretation)  in  the  sense  that  the  theoretical  growth  rates  increase  also,  as  described  in 


the  previous  paragraph.  Thus  the  prediction  here  is  that,  with  cooling,  instability  may  be 
delayed  but  when  it  occurs  it  is  likely  to  be  more  violent  than  in  the  non-cooled,  e.g.  insu¬ 
lated,  case.  This  is  mainly  for  the  viscous-inviscid  first  modes,  but  similar  findings  emerge 
for  the  inviscid  second  modes  also  as  discussed  later. 

“Moderate  cooling”  of  the  surface,  and  its  influence  on  the  compressible  viscous  modes,  is 
dealt  with  in  § §2  and  3  below,  based  on  the  theory  in  Smith  (1987).  The  new  flow  structure 
produced  is  examined  first,  followed  by  its  instability  properties.  The  flow  structure  is  of 
the  compressible  Rayleigh  inviscid  kind  across  the  majority  of  the  boundary  layer,  but  in  a 
quasi-steady  form,  which  admits  pressure-displacement  interaction  with  the  viscous  sublayer 
response.  Unlike  the  non-cooled  TS  instability  modes,  which  in  supersonic  flow  must  be 
directed  outside  the  wave-Mach  cone,  i.e.  be  effectively  subsonic,  the  present  surface-cooled 
ones  can  be  unstable  for  any  obliqueness  angle,  including  the  two-dimensional  (2D)  case. 
This  and  other  factors  of  the  link  between  the  cooled  and  the  non-cooled  TS  modes  are  also 
explored  in  §3,  along  with  the  influence  of  cooling  the  surface  beyond  the  moderate  stage. 
There  the  spatial  growth  rates  become  still  more  accentuated,  for  any  Mach  number,  as  the 
viscous-inviscid  disturbance  is  concentrated  relatively  near  the  surface,  and  flow-structural 
change  re-occurs.  We  devote  a  separate  section,  §4,  to  the  resulting  stage  of  “severe  cooling”, 
to  emphasize  its  novel  nature.  This  concerns  principally  the  intrusion  of  compressibility 
effects  into  the  dynamics  of  the  thin  viscous  sublayer.  Although  the  work  overall  is  mostly 
on  viscous  modes,  pure  inviscid  modes  are  addressed  in  §5  as  a  result  of  the  behavior  found 
earlier  in  the  moderate-cooling  regime.  Once  again,  cooled-surface  modes  which  appear  to 
be  distinct  from  those  found  in  earlier  studies  are  implied,  for  any  Mach  number  A/qo,  with 
the  heat-transfer  sublayer  playing  an  active  part  throughout.  The  new  inviscid  modes,  like 
the  viscous  ones  of  §§2-4,  all  have  relatively  small  phase  speeds,  and  they  include  some  which 
radiate,  i.e.  provoke  outward-traveling  waves  at  the  edge  of  the  boundary  layer. 

Nondimensional  variables  are  used  such  that  the  freestream  velocity  is  ( u ,  v,w)  =  (1,0,  0), 
in  cartesian  coordinates  (x,y,z)  where  the  airfoil  chord  is  unity.  The  nondimensional  pres¬ 
sure,  density,  temperature  and  viscosity  are  written  p,  p,  T,  p  in  turn,  equal  to  px,  1,  1,  1  in 
the  freestream,  where  px  =  1/(7A/^),  and  t  denotes  the  usual  convectively  scaled  time. 
The  basic  boundary  layer  then  has  0(Re~1^2)  characteristic  thickness  where  Re  is  the  global 
Reynolds  number  produced  above  and  the  boundary-layer  properties  are  defined  by  the  pro¬ 
files  (u,  T,  p)  —  (Uo,  T0)  po)(y),  with  y  =  Re~x/2y ,  at  the  local  station  x  =  x0,  z  —  z0  under 
consideration.  In  §§2-4  nonlinear  viscous-inviscid  interactions  are  considered  first  since  the 
understanding  of  flow  transition  is  the  ultimate  theoretical  goal,  although  linear  solution 
properties  are  then  treated,  as  a  first  study,  since  this  type  of  approach  for  surface  cooling 
appears  to  be  a  fresh  one.  The  corresponding  nonlinear  inviscid  versions  in  §5  can  likewise 
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be  set  up.  Further  discussion  is  presented  in  §6,  along  with  some  comparisons  with  the 
experiments  of  Lysenko  and  Maslov  (1984)  and  related  comments. 

2.  THE  MODERATE-COOLING  REGIME 

For  order-one  values  of  the  temperature  ratio  T^/T^(=  Tw )  the  compressible  three- 
dimensional  TS  modes  are  controlled  by  the  triple-deck  structure  as  described  by  Smith 
(1987),  as  long  as  the  other  parameters  involved  are  not  too  extreme.  That  structure  pro¬ 
vides  the  springboard  for  the  present  theory  on  surface  cooling  and  so  its  underlying  form  is 
noted  here.  The  expressions  for  the  flow  solution  in  the  three  decks  are: 


upper  deck,  y  —  Re  3/8K\m  1^2yt  and  ) 

(2.1a) 

[u,  v,  w,p]  =  [1,  0, 0,Poo]  -f  Re  1^4Cl^4X1^2m  1/,4[u^2\  vn}^2T^2\  w^2\  p(2)];  j 
main  deck,  y  =  Re~4^2y  and 

►  (2.16) 

{u,v,w,p}  =  [U^,Q,Voo\  +  0[Re-'l\Re-ll4,Re-'l4,Re-ll4)-  . 
lower  deck,  y  =  Re~5^sm^1^8C5^T^2 \~3/4Y*  and 

[u,v,w,p)  =  {Re-4l*Cxl&T4l2\4l4m~4l8U\Re-y*C3l*T4l2\3/4rnl/*V',  (2.1c) 

i2e-1/8cl/8ri/2A1/4rn-1/8WA.,  Re-lHC'/4\l/2m-l/4 />•]. 
throughout,  the  streamwise  and  spanwise  length  scales  and  the  time  scale  are  given  by 

(x-x0,z-  z0)  =  Re-^K^X^Z*),  t  =  Re-l/4m-l/4TwCl/4\-3/2t\  (2.1  d) 

Here  m  —  \Ml0  —  1|,  and  P  and  —A  are  the  unknown  pressure  and  displacement  variations 
(independent  of  y,Y*),C  is  the  Chapman  constant,  A(x0,zo)  is  the  0(1)  skin-friction  coef¬ 
ficient  of  the  incident  compressible  boundary  layer  and  (xo,Zo)  is  the  local  position  under 
consideration,  while  K\  =  C3^8T^2\~5^  4m  3/,a.  Similar  expansions  hold  for  the  temperature 
T,  the  density  p  and  the  viscosity  p.  As  a  result  the  overall  nonlinear  three-dimensional  TS 
behavior  here  is  governed  by  the  lower-deck  equations  and  boundary  conditions, 


u'x.  + 1/;.  +  wz.  =  o, 


(2.2a) 


U*.  +  U'U'X.  +  VU^.  +  W'U'Z.  =  -P'x.{X\Z\t')  +  Uy*  y  •  i 

w;.  +  u*w*x.  +  vw$.  +  w*wz.  =  -p*z.{x*,z\c)  -\  w^.r., 

u*  =  V*  =  W  =  0  at  Y*  =  0, 

U'  ~  Y*  +  A*{X\  Z\  r),  W*  ->  0,  as  Y*  ->  oo, 


(2.2  b) 

(2.2c) 

(2.2  d) 

(2.2 e) 
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coupled  through  pressure-displacement  interaction  with  the  upper-deck  linearized  potential 
flow  problem  (for  supersonic  and  subsonic  flow  respectively) 

PC  -  l)(a*.  T  o  -  =  0,  (2.2/) 

where 

p^2)  ->  0  as  yu  -»  oo;  p<2)  ->  P*,  ->  A*x.x.  as  yu  ->  0  +  .  (2.2 y) 

Suitable  constraints  at  large  |A"*|,  |Z*[,  are  also  assumed,  and  imposed  forcing  can  be  incor¬ 
porated  in  either  or  both  of  the  parts  (2.2a-e),  (2.2f,g).  Various  nonlinear  features  of  the 
unsteady  viscous-inviscid  system  (2.2)  for  nonlinear  TS  evolution  are  addressed  by  Smith  & 
Walton  (1989),  Smith  (1988,  1989)  among  others,  while  the  linearized  version  corresponding 
to  relatively  small  compressible  TS  waves  is  studied  by  Smith  (1987).  Here  we  follow  the 
same  procedure  of  setting  up  the  nonlinear  stage  first  and  linearizing  afterwards,  in  dealing 
with  surface-cooling  effects. 

The  main  impact  of  surface  cooling,  i.e.  small  Tw,  is  felt  partly  through  the  Tw  factors 
themselves  in  the  scales  above  and  partly  through  the  alterations  in  the  mean-flow  profiles 
Vo,  To  near  the  surface,  specifically  in  the  large  heat  transfer  and  skin  friction  produced  there. 
This  can  be  seen  in  two  ways.  One  is  by  means  of  the  IIowarth-Dorodnitsyn  transformation, 
y’  =  lonely,  since  U0,T0  are  regular  in  terms  of  the  transformed  variable  y'.  So  if  Tw  is 
reduced  then  initially  T0  <x  y'  at  small  y'  with  an  0(1)  constant  of  proportionality.  Hence 
the  transformation  implies  that  y  oc  j0  T0dy'  is  quadratic  in  y',  and  therefore  the  temperature 
profile  T0  is  tending  to  zero  like  the  square-root  of  the  normal  distance  y.  The  same  applies  to 
the  velocity  Vo.  That  implies  large  temperature  and  velocity  gradients  developing  near  the 
surface.  The  precise  size  of  these  is  determined  by  a  sublayer  produced  where  y  is  0(T2),  or 
y'  ~  Tu,,  and  the  low  surface  temperature  becomes  comparable  with  the  linearly  decreasing 
(in  terms  of  y')  contribution  above,  so  that  T0  «  Tw  +  y'  times  0(l),i.e.Tw  +  0(1)  times 
[y  +  0( l)j1/2.  It  follows  that  the  surface  heat  transfer,  dT0/dy  at  y  =  0,  is  0(T~l).  Likewise, 
the  surface  shear  stress  is  0(Tfl).  The  alternative  is  to  consider  the  basic-flow  momentum 
and  heat-flux  equations,  which  are  dominated  by  their  viscous  contributions  near  the  surface 
at  small  y,  i.e.  by  3(/i05V0/ dy) / dy  =  0  and  d(y,0dT0/ dy)/ dy  =  0  respectively.  Since  //0  oc  T0 
for  a  Chapman  fluid  (see  also  later)  the  second  balance  above  integrates  to  give  T0dT0/dy 
uniform  in  y  and  hence  T02  is  linear  in  y;  so  again  the  almost  square-root  response  in  the 
temperature  T0  and  the  large  0(T~l)  heat  transfer  dT0/dy  are  implied  at  small  y  and  Tw. 
Similar  results  for  Un  stem  from  the  momentum  balance  above.  The  effect  of  surface  cooling, 
then,  is  to  induce  a  new  sublayer  in  which,  for  certain  0(1)  constants  C\,c2, 

V  ^  T2J,  To  :  Tw{  1  1  2cly)'/2,  U0  -  Tw^{(  1  +  2 ny)'*2  -  1},  (2.3a) 

<"i 
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and  in  consequence  the  scaled  heat  transfer  and  skin  friction 

(dTo/dy)y=o  ~  T~l ci,  A  =  (dUo/dy)^  ~  T?c2i  (2.3 b) 

are  increased  substantially. 

The  implications  for  the  TS  waves  follow  directly  from  the  increase  in  A  in  (2.3b)  and  the 
decrease  in  Tw.  Thus  the  TS  length  scales  of  (2.1a,c,d)  all  shrink  now,  as  might  be  expected 
physically,  and  in  particular  the  streamwise  length  and  the  maximum  normal  distance  both 
decrease  like  T^2 .  Similar  deformations  occur  in  the  typical  lower-deck  thickness  and  in 
all  the  velocity  and  pressure  scales,  for  instance,  and  eventually  a  new  flow  structure  must 
come  into  play  as  Tw  continues  to  drop.  Various  options  present  themselves  but  the  first 
distinct  new  structure  is  found  to  arise  when  Tw  is  decreased  to  the  order  Re~l/U ,  at  which 
stage  the  upper-deck  extent  becomes  coincident  with  the  main  boundary  layer’s,  as  does 
the  streamwise  length  scale.  So  novel  elliptic  effects  must  then  enter  the  boundary-layer 
dynamics. 

In  the  new  stage  of  “moderate  cooling”,  therefore,  we  have 

Tw  =  Re~l/12TW,  (2.4) 

say,  with  Tw  typically  of  order  unity,  and  a  three-tier  structure  comes  into  operation.  See 
Figure  1.  The  scales  and  balances  involved  are  suggested  by  the  ones  mentioned  earlier  but 
combined  with  (2.4).  Across  most  of  the  boundary  layer,  the  main  tier  where  y  is  0(1)  is 
given  by 


[u,v,w,p,pfT  j  =  [Uo,  9,  OjPoo,  po,  Tq]  +  Ryll4pA2cl[u,v,w,p,p,T} 


nth  c3  —  (C/Ap2)1Oj  while  in  the  buffer  tier  where  y  —  Re  2^2y 


[u,v,w,p,p,T]  =  [Re  l!uu  +  Re  1/6c3iu ,  -Re  1/3c3Axu,Re  l^X2c23D /u, 

Poo  +  Re-1/4p\2c23P,Rel/l2p,Re-l'l2fw}, 
to  leading  order,  with  Dx  =  -Pz,  and  in  the  lower  tier  y  =  Re  3'/'1c3 Y  with 


[u,  v,  w 


,p,p,T }  =  {Re-l/*\c3U,  Re-*/l2\clV,  Re^/6\c3W, 


Poo  +  Re-llAp\2c\P,  Rel/np}  Re-l'12fjr  (2.7) 

Here  the  extra  factors  0,  A(=  c2/Tw),  p(~  Twl)  present  are  all  0(1)  and  are  included  to  pro¬ 
duce  a  normalized  equation  set  below,  taking  account  of  the  scaled  surface  density,  viscosity 
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and  shear  values  which  would  otherwise  appear  in  the  lower-tier  system  for  example.  The 
streamwise  and  spanwise  length  and  the  time  scales  now  are  defined  by 

Re~l/3t 

(x  -x0,z-  z0)  =  Re-^2(X,  Z),  t  =  tZL - .  (2.8) 

Ac3 

Hence  the  main-tier  properties  are  controlled  by  the  linearized  quasi-steady  inviscid  balances 


Poux  +  PxUq  +  poVy  +  pot  v  +  Pou>z  =  0,  (2.9 a) 

Uqux  +  v  U0y  =  -px/po,  (2.9  6) 

Uo(vx,wx)  =  ~(Py,Pz)/po,  (2.9c) 

Uo(poPx  -  'ypoopx )  =  IPoopOyV,  (2.9 d) 

leading  to  the  compressible  pressure  equation 

Pyy  +  Pzz  +  (^  —  M2)px  x  =  2Mq -Py!  Mo,  (2.10a) 


with  effectively  zero  phase  speed.  Here  M0  =  U0M~^  jT^2 .  The  constraints  on  p  are 


p  — »  0  [or  outgoing  waves]  as  y  — >  oo,  (2.106) 

2 

P  ~  P  +  -rAxxy212  as  y  ->  0+,  (2.10c) 

together  with  appropriate  boundedness  conditions  far  upstream  and  downstream  (see  below), 
where  (2.10b)  reflects  the  match  with  the  oncoming  undisturbed  free-stream  motion  and 
(2.10c)  is  necessary  for  merging  with  the  buffer  tier.  The  3/2-dependence  in  (2.10c)  comes 
from  analysis  of  (2.10a),  given  that  the  basic-flow  profiles  Uq,Tq  exhibit  the  square-root 
dependence 

Uo  -  c2(2/c1)I/¥/2,  T0  ~  (2Cl)1/2y1/2,  as  y 0  +  ,  (2.11) 


consistent  with  (2.3).  The  0(1)  parameter  r  introduced  in  (2.10c)  is  defined  by 

c2/321/2r2 


T  = 


3/2 


<71/3 


(2.12a) 


for  consistency  with  the  buffer  tier  below.  Again,  the  associated  velocity  perturbations 
respond  in  the  form 

(u,v,w)  ~  (y~1/2,y1/2,0(  1)),  (2.126) 


at  small  y,  agreeing  both  with  (2.10c)  and  (2.6).  The  buffer  tier  of  (2.6)  is  also  governed  by 
predominantly  inviscid  linear  dynamics  but  it  accommodates  the  basic  inner-layer  behavior 
(2.3a),  thereby  smoothing  out  the  3/2-dependence,  since  the  velocity  profile 

fi(»)=-[(i2  +  2c,j)1'J-7U  (2.13) 

Cl 
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is  regular.  The  final,  lower,  viscous  tier  of  (2.7)  then  yields  the  three-dimensional  unsteady 
nonlinear  boundary-layer  equations, 


U  x  +  Vy  +  Wz  —  0,  (2.14a) 

Ui+UUx  +  VUY  +  WUz  =  -Px(X,Z,i)  +  UYY,  (2.146) 

Wi  +  UWx  +  VWY  +  WWz  =  -Pz(X,  Z,  t)  +  WYY,  (2.14c) 

subject  to  the  conditions 

u  =  V  =  W  =  0  atY  =  0,  (2.Ud) 

U  ~  Y  +  A(X,  Z,  t),  W  —+  0,  as  Y  — >  oo,  (2.14e) 


for  no  slip  at  the  surface  and  to  match  with  the  buffer  tier’s  response  in  turn. 

The  nonlinear  system  describing  moderate  cooling  is  therefore  given  by  the  viscous- 
inviscid  interaction  of  (2.10a-c)  combined  with  (2.14a-e).  A  match  may  be  verified  with  the 
earlier  structure  of  (2.1)  at  large  r,  corresponding  to  large  Tw,  and  certain  extra  aspects  of 
this  are  examined  in  the  next  section.  Of  most  concern  here  are  the  properties  of  the  new 
cooled-surface  TS  system  at  0(1)  values  of  r  and  the  effects  of  the  inviscid  but  rotational 
elliptic  part  (2.10)  of  the  interaction,  which  represents  the  major  new  feature  present  in  the 
moderate-cooling  regime.  Further,  small-disturbance  properties  can  be  of  much  interest,  for 
practical  and  other  reasons  including  the  fact  that  the  length  scales  and  hence  the  spatial 
growth  rates  expected  are  now  comparable  with  those  of  pure  inviscid  modes,  and  so  we 
address  here  the  linearized  instability  version  in  which  (U,  V.  W.  Py  A ,  p)  are  taken  to  be  small 
perturbations  of  the  basic  undisturbed  state  (Vj  0, 0, 0,  0, 0).  This  has  the  effect  of  keeping 
(2.10a-c)  intact  but  with  the  lower-tier  response  becoming  linearized.  In  consequence  cool- 
surface  TS  waves  oc  exp(i(aX  +  (3Z  —  Q  t ))  are  possible  with  wavenumbers  a,/3  and  frequency 
Q,  in  which  case  we  are  left  with  solving  the  quasi-steady  compressible  Rayleigh  equation 


subject  to 


^p'-[aJ(l-A/o!)  +  ^lp  =  0, 

(2.15 a) 

►  0  [or  outgoing  waves]  as  y  — >  oo, 

(2.156) 

~  P  -  ^r^Ay3/2  as  y  ->  0+, 

(2.15c) 

y,  (6  =  -i1/3f2/a2/3  ,  k  =  Ai(q)dqSj  . 

(2.1 5d) 

and 

P  (ra)-1/3 

A  ~  (1  +  /?2/a2) 

In  (2.15a)  the  prime  denotes  a  y  -  derivative,  while  (2.1 5d )  stems  from  the  viscous  linear 
lower-tier  solution  in  terms  of  the  Airy  function  Ai.  The  properties  of  the  moderate  -  cooling 
system  (2.15)  are  investigated  in  the  next  section. 


7 


3.  INSTABILITY  PROPERTIES  FOR  MODERATE  COOLING 


Solutions  of  the  moderate  -  cooling  problem  (2.15)  were  obtained  (a)  by  computation,  (b) 
analytically  for  small  t,  (c)  analytically  for  large  r,  connected  with  relatively  high  frequencies, 
and  these  are  described  in  turn  below.  In  practice,  some  of  the  features  found  in  (b),  (c) 
were  used  to  guide  the  computations  in  (a),  and  vice-versa. 

3(a).  Computational  Study 


For  moderate  values  of  the  Mach  number  we  chose  the  mean-flow  profiles  U0  and  T0  to 
have  the  following  representative  forms: 


T0  =  1  -  e“"1/2, 


(3.1a) 

(3.16) 


i.e.  they  have  the  required  square  -  root  dependencies  as  y  — >  0  +  ;  the  decay  as  y  — >  co  is 
of  relatively  little  consequence,  we  note.  Then  M0  =  UoM^/T^2  =  M^l  —  and 

the  system  (2.15)  becomes 


1  p~  ylp 

?"  -  7jT'n- - wtV  -  (a2(l  -  ML (1  -  e-«'n)  +  P2)p  =  0, 

^  1  —  e  y 


with  the  boundary  conditions 


and 


p  — >  0  (or  outgoing  waves)  as  y  — ►  oo, 

2 

p  ~  P - r^Ay3/2  as  y  0  +  , 

3 

P  (za)"1/3  Ai'((o) 


(3.2a) 

(3.26) 
(3.2 c) 

(3.2  d) 


A  (l+/?2/«2)  «(fo)  ’ 
where  £0  and  k  are  defined  by  (2.15d).  For  y  <<  1  we  may  obtain  two  independent  solutions 
of  (3.2a),  with  multiplicative  constants  aj,  a2, 


2  Ml 


P  =  «,<1  +  (P2  +  c2)f  -  \(p2  +  a>(l  +  2Ml))t^  +  (A(^  +  a1)  +  +  •  •  ■} 


+«2{y3,!  -  +  l^t'2  -  +  ••■}.  for  y  «  i. 


(3.3) 


One  numerical  method  used  to  solve  the  eigenvalue  problem  (3.2)  is  to  start  the  integration 
procedure  at  a  small  value  of  y  =  6,  say.  Then  for  non-neutral  solutions,  for  example,  we 
can  fix  ft  and  guess  a,  complex,  and  find  two  independent  solutions  for  p  corresponding 


8 


to  (01,^2)  —  (1,0),  (0,1)  using  a  standard  4th  order,  4  stage  Runge-Kutta  scheme.  The 
integration  is  continued  to  a  large  value  of  y,  say,  where  p  ~  exp[-(a2(l  -  AO  +  /52)1/2y) 
and  the  two  solutions  pl  and  p2  are  combined  in  the  form 

2 

P  =  PPi  -  ^ra2p2,  (3.4) 

where  P  is  defined  by  (3. 2d).  Without  loss  of  generality  we  may  choose  A  —  1.  This  system 
is  iterated  until  the  condition 

p’(yoo)  +  ^p(yoo)  =  0  (3.5) 

is  satisfied,  where  A  =  (a2(l  —  A/£,)  +  /52)1/2  with  Re{ A)  >  0.  In  another  method,  to  avoid  the 
possibility  of  exponential  growth  in  the  numerical  solution  on  the  interval  (6,  y<x>),  the  solution 
of  (3.2a)  can  be  found  first  from  y  =  8  to  y  —  y0  (y0  =  1,  for  instance),  choosing  .4  =  1  and 
using  the  asymptotic  behavior  of  (3.3)  as  initial  conditions.  Additionally  a  solution  of  (3.2a) 
may  be  found  from  y  =  y^  to  y  =  yo  using  the  exponential  decaying  behavior  at  y  =  y00) 
and  choosing  the  constant  multiplying  the  solution  so  that  the  two  solutions  are  continuous 
at  y  =  y0.  The  condition  that  p  and  p  are  continuous  at  y  =  yQ  is  sufficient  to  determine 
the  complex  eigenvalue  a  for  fixed  /5,  ¥«,,  r  and  Pl.  In  some  instances,  to  aid  convergence, 
a  further  method  was  used  where  it  was  found  useful  to  map  (0,oo)  onto  (0,  1)  by  use  of 
the  transformation  y  —  tan(7rs/2).  The  method  of  solution  in  this  case  can  be  either  of  the 
two  methods  described  above,  since  the  asymptotic  behavior  of  the  transformed  equation  of 
(3.2a)  as  s  — >  0T  and  3  — >  00  is  unchanged  (simply  replace  y  by  tan(7rs/2)).  Each  method 
described  above  was  found  to  give  results  consistent  with  the  other  two.  Each  scheme  was 
used  with  different  values  of  6  and  the  step  length  for  the  Runge-Kutta  integration  scheme 
to  test  for  convergence. 

Neutral  and  unstable  solutions  of  (3.2)  were  found  by  using  the  numerical  methods  de¬ 
scribed  above  and  are  shown  mainly  in  Figure  2;  see  also  in  Figures  3,4  later.  The  behavior 
of  the  solutions  was  followed  as  the  (real)  frequency  Pl  was  increased.  The  neutral  solutions 
were  obtained  for  two-dimensional  and  three-dimensional  disturbances.  It  was  found  useful 
to  plot  the  neutral  curves  for  fixed  values  of  (3 /a.  Figure  2(a)  shows  the  neutral  values  of  r 
and  Pl  for  Mx  =  0.5  and  (3/ a  =  0, 1,2,  4.  We  see  that  r  decreases  as  Pl  increases.  The  neutral 
solutions  for  supersonic  flows  difur  depending  on  whether  /5/a  —  ( M ^  —  1)5  is  positive  or 
negative.  (The  situation  with  /5/a  =  ( M —  1)j  is  not  addressed  here.)  Figure  2(b)  shows 
neutral  values  of  r  and  Pl  for  =  2  and  /5/a  =  1.75,  1.5,  1.25, 1,0.  For  /5/a  >  (M^,  —  1)? 
the  solution  for  p  is  real  and  decays  exponentially  as  y  — >  00.  For  /5/a  <c  we 

allow  outgoing  traveling  waves  and  so  p  is  complex  and  has  a  purely  oscillatory  behavior  as 
y  —>  00.  These  two  different  types  of  solution  produce  different  forms  of  the  neutral  curves 
as  shown  in  Figure  2(b).  For  /5/a  >  ( —  1)?  (/5/a  =  1.75)  there  is  no  solution  for  large 
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values  of  both  r  and  Q.  For  p/a  <  (M^  —  1)5  there  exists  a  cut-off  value  of  Q  for  each 
value  of  P/a  below  which  there  are  no  solutions.  As  Q  increases  from  the  cut-off  value, 
on  the  upper  curve,  r  reaches  a  local  maximum  and  then  decreases  to  a  local  minimum 
as  fi  increases  further.  Beyond  the  local  minimum  the  curve  has  no  more  turning  points 
and  r  continues  to  increase  as  increases.  We  note  that  the  upper  curve  for  the  2D  case 
P/a  =  0  does  have  this  behavior  but  because  of  the  scale  of  Figure  2(b)  it  is  not  immediately 
apparent.  The  lower  curves  become  almost  coincident  as  Q  increases,  including  the  cases 
P/a  >  (A/£,  —  I)?  for  which  r  continues  to  grow  as  Q  becomes  increasingly  small.  Thus  we 
see  that  there  exist  neutral  (and  hence  unstable)  solutions  of  (3.2),  over  a  range  of  r,  for  all 
values  of  fi  for  p/a  >  (A/£,  —  1)5  or  for  any  subsonic  case,  but  there  is  a  cut-off  value  of  Q 
for  each  value  of  P/a  for  which  P/a  <  (A/£,  —  1)5  where  no  neutral  (or  unstable)  solutions 
exist  for  il  below  this  value.  Figure  2(c)  shows  neutral  values  of  a  and  Q  for  M ^  =  2  and 
Pj a  —  1.5,  1.25,  1, 0.  The  curves  for  P/a  >  —  1)5  for  the  supers  onic  case  and  for  any 

subsonic  case  are  also  shown.  They  are  only  present  on  the  upper  curve  wrhere  the  different 
solutions  are  indistinguishable.  For  P/a  >  -  1)5  and  also  for  <  1  the  pressure 

p  is  real,  hence,  from  (3. 2d)  P  must  be  real  also.  This  implies  the  well  known  result  that 
£o  =  — 2.298r5  and  i~  >  Ai' ((0) / k((o)  =  1-001.  The  definition  of  £0  from  (2.15d)  determines 
a  for  these  cases.  We  find  that  a  =  (0/2.298)5.  Hence,  for  the  cases  where  this  solution  is 
valid,  a  is  independent  of  Moo,  P  and  r.  The  two  branches  of  the  neutral  curves  shown  in 
Figure  2(c)  correspond  to  the  asymptotic  solutions  described  later  in  sections  3(b)  and  3(c). 

The  effects  of  increasing  the  cooling  of  the  surface  were  investigated  by  decreasing  r  for 
two-dimensional,  three-dimensional,  subsonic  and  supersonic  (for  moderate  values  of  Moo ) 
disturbances.  The  trends  were  found  to  be  the  same  in  all  of  these  cases,  namely  that 
decreasing  r  increases  the  growth  rate  of  the  unstable  solution.  This  is  shown,  for  example, 
in  Figure  2(d)  for  =  0.5,  P  =  1  and  r  =  0.1,0.05,0.02,  where  the  growth  rate  -7m(a) 
is  plotted  against  fh  We  observe  that  the  growth  rate  increases  significantly  as  r  decreases 
(i.e.  when  the  wall  is  cooled  more).  This  suggests  a  large  destabilization  of  the  basic  flow 
as  the  wall  cooling  is  increased.  The  results  for  A/oo  —  0.5  and  p  =  0  for  r  =  0.1,0.05,0.02 
are  graphically  identical  to  the  respective  results  for  P  —  1  shown  in  Figure  2(d).  The  real 
part  of  a  is  very  large  for  these  solutions  for  small  values  of  r.  Figure  2(e)  shows  the  growth 
rate  —  Irn(a)  as  a  function  of  fi  for  Mx  —  2,  /?  =  1  and  r  =  0.1,0.05,0.02.  The  solutions 
have  the  same  behavior  as  those  shown  in  Figure  2(d)  for  Moo  —  0.5  with  the  growth  rates 
appearing  to  be  slightly  smaller  for  M oo  —  2  compared  to  those  for  A /^  =  0.5.  Again,  the 
corresponding  solutions  for  A/CXJ  2  and  p  -  0  are  virtually  the  same  as  those  for  A/^  —  2 
and  P  1  for  the  values  of  r  shown  in  Figure  2(e). 

For  large  values  of  r  the  solutions  approached  those  of  Smith  (1987)  as  expected  [.-.ee 
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also  (c)  below].  For  subsonic  flows  as  t  increases  the  solutions  for  Im(a)  have  tne  same 
behavio’-  as  those  shown  ;n  Figure  2(d)  with  the  growth  rates  continuing  to  decrease  as  r 
increases.  This  behavior  is  found  to  hold  for  very  large  values  of  r.  This  is  illustrated  in 
Figure  2(f)  where  —Im{a)  is  plotted  against  f2  for  A/^  —  0.5,  [3  —  0  and  r  =  0.4,  0.5,1. 
If  the  corresponding  curves  for  f3  =  1  were  to  be  added  to  Figure  2(f)  they  would  be 
indistinguishable.  Figure  2(g)  shows  —  /777(a)  as  a  function  of  Q  for  Mx  —  0  5  and  j3  -  0 
and  (3  =  1  for  r  —  10,20.  The  solutions  for  (3  ~  \  and  (3  =  0  are  now  distinct,  with  the 
solutions  for  (3  —  0  having  larger  growth  rates,  the  difference  increasing  as  r  increases.  We 
note  that  as  fi  increases  from  zero  the  three-dimensional  mode  occurs  first. 

For  supersonic  flows,  as  expected  from  the  neutr  d  curves  shown  in  Figure  2(b),  the 
structure  of  the  eigenvalue  a  changes  as  r  increases.  An  example  of  this  is  shown  in  Figure 
2(h)  where  —  Im(a)  is  plotted  against  for  Moo  —  2  and  (3  =  0  for  r  =  0.4,  0.5,  1.  For 
0  <  r  <  0.4  there  is  one  neutral  point  at  a  small  value  of  Q.  For  0.4  <  r  <  1  there  are  3 
neutral  points  and  for  r  >  1  only  the  neutral  point  at  the  largest  value  of  Q  persists.  Figure 
2(i)  shows  — Im(a )  as  a  function  of  for  larger  values  of  r,  namely  for  Moo  =  2,  (3  ~  0  and 
r  —  5,  10,  17.  For  the  range  of  r  shown  in  Figures  2(h)  and  2(i)  the  solutions  for  (3  —  1  are 
very  close  to  those  shown  for  /?  =  0.  As  r  increases  the  growth  rate  continues  to  decline. 

To  summarize  the  results  of  the  computational  study  for  moderate  values  of  M ^  we 
can  make  the  following  general  statements.  As  the  amount  of  wall  cooling  is  increased 
(he.  as  t  is  decreased)  the  growth  rates  increase  for  two-dimensional  and  three-dimensional 
perturbations,  irrespective  of  whether  the  flow  is  subsonic  or  supersonic.  Moreover,  as  r 
is  increased,  for  subsonic  flows  no  significant  three-dimensional  effects  are  noticeable  until 
r  >  10,  with  the  growth  rates  here  being  larger  for  two-dimensional  perturbations.  On  the 
other  hand,  for  supersonic  flows  with  (3 / a  <  there  is  very  little  difference  between 

the  solutions  for  two-dimensional  and  three-dimensional  flows.  Again,  for  a  fixed  value  of  r 
the  growth  rates  for  subsonic  flow  are  larger  than  those  for  supersonic  flow. 

We  also  note  that,  here  and  below,  multiple  modes  may  be  present  but  these  have  not 
been  detected  as  yet. 

For  large  values  of  A/^  the  mean  flow  becomes  different  in  structure.  In  particular,  from 
the  IIowarth-Dorodnitsyn  transformation,  T0  ~  M^  so  that  y  ~  A/£,,  while  Uo  ~~  l^(l)-  We 
do  not  present  solutions  for  large  Mach  number  but  below  is  an  outline  of  a  representative 
problem  to  be  solved  in  this  case.  For  large  M^  we  write 

y  VAfl,  (3.6) 

«]  mV),  A/^r0(F),  A7„(F),  Mja], 

1 1 


and  then 


|F),  I  i i,  A  /o , 


(3.7) 


instead  of  (3.1a,b),  with  the  over-barred  variables  now  being  typically  of  order  one.  The 
equation  (2.15)  therefore  becomes 

PTY  -  -  Ml)  +  p‘\p  0.  (3.8) 

for  0  <  Y  <  Y oo,  where  the  scaled  boundary-layer  edge  Y is  finite  in  this  regime,  M0  = 

_  _ _ 1  ,  _ 

UqT0  5,  and  /3  =  M^2/?.  The  boundary  conditions  near  the  cooled  surface  likewise  remain 
intact,  to  leading  order,  with  z ,  P ,  A  also  scaled  appropriately  with  respect  to  M <*,,  taking 
account  also  of  the  extra  factors  such  as  C\ ,  c2  involved  in  the  scales  of  section  2.  In  particular, 
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the  typical  wall-cooling  factor  r  required  now  increases  like  M .  A  logarithmically  thin  extra 
layer  near  Y  =  Y adjusts  the  solution  from  T0  being  O(M^)  (as  above)  to  being  unity  as 
required  in  the  free  stream,  and  so  the  profiles  U o,  T0  are  to  satisfy  the  conditions  C/ 0  — »■  1 , 
T0  -»  0  as  F  ->  Foo-  and  U0  oc  F’,  T0  oc  Y*  as  F -*  0+. 


(b)  Small  r:  towards  severe  cooling 

In  line  with  the  computational  results,  it  is  suggested  that  for  small  values  of  r  the 
disturbance  y  -  scale  contracts.  So  U0,T0  are  approximately  like  y1^2,M0  is  small  of  order 
y1/*  for  Aloo  of  0(1),  cf.  later,  and  therefore  (2.15a)  suggests  that  a  is  0(j/r~1),  as  is  (3 . 
The  unknown  y-scale  is  then  fixed  in  terms  of  r  by  comparing  the  typical  size  of  F/A  from 
(2.15c)  with  that  inferred  from  (2.15d)  ,  given  that  £0  is  expected  to  remain  typically  0(1) 
for  instability.  The  implied  small-r  response  therefore  has 

p  =  p(pH — »  y  =  r6/5y,  (3.9a) 

(Q,  £0  —  r~6/5(d,/3)  +  •  •  • ,  n  =  --4/50  +  • ,  (3.96) 

P  =  p  +  --  ,  A  -  r  2/5  i  1  •••.  (3.9c) 

This  reduces  (2.15)  to  the  equation 

hi  ~  ft  -  (AJ  +  hi>  =  0,  (3.10a) 

for  p,  subject  to 

p(oo)  =  0,  p  -  P  —  -a2Ay3//2  as  y  — >  Of,  (3.106,  c) 

3 

and 

PI  A  -  (1  i  f)2/a2)l{iA)  1/3Ai((0)/k{( o),  (3  10d) 

where  ( 0  =  q:!//3.  Here  the  solution  of  (3.10a,b)  can  be  obtained  analytically, 

p  --  til y3//4 Ks(y)  with  y  =  (q2  1/5 2)1/2y,  (3.11a) 
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d\  is  a  constant,  A"3/4  is  the  Bessel  function  in  the  standard  notation,  and  then  (3.10c) 
requires  that 


P/A  = 


27TC*2 

(r(?))2(a2  +  /32)3/4 


(3.116) 


from  the  behavior  of  the  Bessel  function  at  small  y.  So  the  coupling  with  (3 . 1  Od )  leads  to 
the  reduced  eigenrelation 


Ai'Uo) 

*(io) 


2ft 

(f (W 


(3.11c) 


between  a,  Cl.  The  solution  for  Re(a)  and  7m( ot)  from  (3.11c)  are  plotted  in  Figures  3(a) 
and  3(b),  respectively,  for  (3—1  and  r  =  0.1,0.05,0.02.  These  are  seen  to  agree  well  with 
the  full  computations  at  small  values  of  r;  see  also  in  Figures  2(d)  and  2(e).  In  order  to 
illustrate  this  agreement  Figures  3(c)  and  3(d)  show  Re(a )  and  Im(a)  for  r  =  0.05,  (3  =  0, 
for  AIqo  =  0.5  and  M =  2  from  the  full  computations  along  with  the  result  of  (3.11c).  The 
agreement  is  seen  to  be  very  close  when  fi  is  not  very  large. 

It  is  interesting  that  the  small-r  response  here  is  independent  of  the  Mach  number,  as 
might  be  expected  for  disturbances  focussed  near  the  surface.  This  is  because  the  effective 
Mach  number  M0  is  small,  so  that  the  contribution  Mq  in  (2.15a)  is  negligible,  although  the 
instabilities  that  are  found  to  exist  do  remain  temperature-  and  velocity-driven  through  the 
ratio  Moy/Mo  operating  in  (2.15a).  The  property  that  these  extra-cooled-surface  instabilities 
exist  at  any  Mach  number,  and  for  both  2D  and  3D  modes,  contrasts  with  the  situation  for 
the  less  cooled  surfaces  corresponding  to  large  r  or  to  the  earlier  regime  of  (2.1),  where,  in 
the  supersonic  range,  instabilities  occur  only  for  sufficiently  oblique  modes  directed  outside 
the  wave-Mach-cone:  see  (c)  below.  Thus  the  cooling  now  destabilizes  a  range  of  previously 
stable  viscous  modes.  Moreover,  the  de-stabilization  is  very  marked,  for  all  the  viscous 
modes,  in  the  sense  that  the  spatial  growth  rates  produced  are  much  larger  than  encountered 
before.  Indeed,  the  growth  rates  appear  to  exceed  those  of  purely  inviscid  modes  (see  also 
section  5).  The  small-r  response  leads  on  to  the  new  structure  investigated  in  section  4. 


(c)  Large  r  and  the  cut-off 

The  clue  to  the  behavior  at  large  r,  associated  with  less  cooled  surfaces,  is  found  in 
high-frequency  properties.  Thus  when  fl  is  large,  at  0(1)  values  of  r  first,  the  outer  solution 
of  (2.15a)  takes  the  form 


[p,  P,  A,  a,  (3,  M0]  [p,  P,  D2/3d,  02/3/?,  0(Cl-1/6yl/4)]  +  ---,y  =  2/3y,  (3.12a) 

concentrated  again  near  the  wall.  So  p  satisfies  the  equivalent  of  (3.10a),  subject  again  to 


(3.10b.c).  Here  however  the  sublayer  relation  is  simply 


P/A  =  a-x(l  +/32/d2)-1, 


(3.126) 


as  )£o|  is  large  for  these  high  frequencies.  Hence  we  have  the  result 

«  =  (r(l)),/3/(2^)2'3. 


(3.12c) 


for  the  2D  example  [the  3D  version  has  an  extra  factor  (1  -f  02 /a2)-1/6  on  the  right-hand 
side],  giving  neutral  stability  at  leading  order.  This  result  is  essentially  an  inviscid  one  since 
the  asymptote  — £0  for  Ai  / k  in  (3 . 1  Od) ,  used  in  (3.12b),  corresponds  to  an  inviscid  limit, 
and,  like  that  in  (b)  above,  the  present  result  applies  for  all  Mach  numbers. 

Closer  inspection  then  shows  that  the  high-frequency  structure  and  results  (3.12b,c) 
continue  to  hold  for  all  r,  provided  that  fi  >>  r  4 A  and  D  >>  r,  the  former  to  keep  |£0| 
large  and  the  latter  to  keep  a  large  (and  y  small),  as  supposed  above.  The  former  restriction 
is  associated  more  with  decreasing  t  since  (b)  comes  into  play  exactly  when  fi  reaches  the 
order  r~4//5  and  thus  provides  a  match.  The  other  restriction  is  connected  with  the  large-r 
behavior  below.  Also,  we  note  that  higher-order  contributions  in  inverse  powers  of  D  bring 
in  the  second  term  in  Ai  / /c,  which  is  a  viscous  and  de-stabilizing  effect,  such  that  there 
exists  instability  for  any  value  of  r,  at  sufficiently  high  frequencies. 

To  complete  the  picture  at  large  r,  then,  we  move  on  to  consider  0  values  of  order  r 
(Figure  4).  The  properties  are  mainly  inviscid  again  since  |£0|  is  large,  in  fact  comparable 
with  Q,  as 

r  ~  Qt  with  a,0  =  0(1),  y  =  0(1),  (3.13a) 

to  leading  order.  The  dominant  equat:  >n  to  solve  is  therefore  exactly  (2.15a),  with  the 
constraint  (2.15b)  still,  but  (2.15c)  becomes 


P  ~  P[1 


(3.136) 


because  the  sublayer  response  gives  here  f }A  —  (a2  +  02)P/a.  The  computed  solution  of 
(2.15a,b)  with  (3.13b)  is  presented  in  Figure  4(a)  along  with  comparisons  with  the  compu¬ 
tations  of  section  3(a)  and  the  solution  from  (3.12).  At  small  f,  corresponding  to  reduced  r 
and/or  larger  D  values,  there  is  a  match  with  (3.12)  as  expected,  for  any  Mach  number.  Fig¬ 
ures  4(b,c,d)  show  solutions  from  the  computations  of  section  3(a)  for  large  r  with  =  2 
and  0  =  0.  In  Figures  4(e,f,g)  we  show  f  as  a  function  of  along  the  neutral  curves  for 
Mao  =  2  and  0  —  0,  l,  1.5,  as  shown  in  Figure  2(b).  We  see  that  for  the  branch  r  >>  1 
and  fi  »  1  f  decreases  as  Q  increases  so  there  will  be  a  cut-off  value  of  f  above  which  no 
unstable  solutions  exist.  As  r  increases  the  solution  becomes  Mach-number  dependent,  and 
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this  is  especially  so  at  large  f  where  the  subsonic  realm  remains  always  unstable  but  the 
supersonic  realm  becomes  stable  within  the  wave-Mach-cone  (3  <  —  l)17"2,  including 

the  2D  case.  The  cut-off  here,  for  the  less  oblique  supersonic  modes,  is  profile-dependent 
but  occurs  at  an  0(1)  crossover  or  neutral  value  of  f.  That  yields  the  predictions  (in  line 
with  the  computations  just  mentioned  and  with  those  of  Figure  2) 


Q  —  neutral  oc  r 
f]  —  max  oc  t 


(3.14a, b) 


for  Moo  >  1  and  waves  directed  inside  the  wave-Mach-cone.  This  behavior  is  illustrated  in 
Figure  4(h)  where  Q  —  max  is  plotted  as  a  function  of  r  for  /3  =  0  and  M x  —  0.5.  The 
frequency  Q  —  0  —  max  is  that  at  which  the  maximum  growth  rate  occurs.  In  contrast, 
analysis  for  the  more  oblique  supersonic  modes  and  for  the  subsonic  ones  shows  that  at  large 
r  the  outer  solution  develops  a  double  structure,  leading  to  the  typical  response 

(3.15) 


for  the  2D  subsonic  case  for  instance.  This  confirms  the  suggestion  of  the  computations, 
that  the  result  (3.14b)  in  fact  holds  for  all  the  supersonic  and  subsonic  modes,  whereas  the 
subsonic  and  more  oblique  supersonic  ones  retain  instability  at  large  f  (i.e.  fi  <<  t),  as 
opposed  to  (3.14a). 

The  double  structure  just  mentioned  leads  to  the  final  phase  for  the  latter  modes,  namely 
when  Q  falls  to  Then  |£o|  decreases  to  0(1)  since 

(a,/3)  ~  r_30(a,,/3.)  with  Q  ~  r_1/2n,,  (3.16) 


from  (3.15),  with  a,,  /?*,  Q,  of  order  unity,  so  the  fully  viscous  sublayer  relation  (2. 1 5d )  is  re¬ 
installed.  Moreover,  the  outer  double  structure  produces  the  two  y-scales  0(r3/4)  and  0(1), 
the  former  giving  a  potential-flow  character  and  the  latter  a  profile-dependent  displacement 
effect,  leading  to  the  relation  P  =  a,Ar~3^  (1  —  for  subsonic  2D  modes  for  instance. 

Coupling  this  with  (2.15d),  then,  we  obtain  the  eigenrelation 


«(&>) 
At  (to) 


(ia.y/3 


(1-Mj)1/2 

a. 


(3.17) 


between  a,,Qt,  again  for  the  2D  subsonic  case  as  an  example.  There  are  two  main  points 
here.  First,  the  relation  (3.17),  which  agrees  with  (3.15)  at  large  Q,  as  required,  produces  a 
maximum  growth  rate  at  an  0(1)  value  (approx.  4)  of  0,  and  neutral  stability  at  a  smaller 
0(1)  value  (0.  approx.  2.3),  thus  giving  rise  to  the  predictions 


Q  —  neutral 
D  -  max 


t-1/2 
T  '  /2 


as  r  — »  oo. 


(3. 18a,  5) 


These  again  appear  to  be  consistent  with  the  computations  at  large  t,  as  shown  in  Figure  4(i) 
where  Vl-max  is  shown  as  a  function  of  r  for  /3  =  0  and  M0 0  =  2.  Second,  the  entire  structure 
in  this  regime,  and  the  results  (3.16)  -  (3.18),  match  directly  with  the  earlier  triple-deck  case 
(see  (2.1)).  For  instance,  the  wavelength  0(i2e_1^2a_1 )  increases  to  0(Re~3^8)  when  Tw  is 
raised  to  0(1),  since  then  Tw  oc  r1^2  — >  0 (Re1//12)  and  a  — >  0(Re~R8)  from  (3.16).  Similar 
matching  holds  with  regard  to  (3.17),  (3.18),  for  2D  and  3D  subsonic  or  supersonic  modes 
(Smith  1987). 

The  splitting  of  the  neutral  curves  for  the  different  cases,  and  the  production  of  two  local 
maxima  (and  mounds  of  instability)  in  certain  cases,  via  (3.14),  (3.18),  reflects  the  emergence 
of  two  major  types  of  instability  when  the  surface  temperature  increases.  These  are  the 
viscous-inviscid  triple-deck  type  (of  (2.1)  ff),  with  which  (3.18)  connects  [see  comments 
above],  and  the  inviscid  compressible-Rayleigh  type,  which  are  implied  by  (3.14).  In  the 
latter  case,  raising  Tw  back  to  0(1)  corresponds  to  Tw  increasing  to  0(/£e1,/12),  hence  r 
to  0(Re O6),  and  so  the  fFs  in  (3.14)  become  0(i?e1/,e)  also,  meaning  that  the  time  scale 
t  becomes  0(Re~1^2),  from  (2.8).  This  t-scale,  the  x-,y-scales  of  order  Re~1^2,  and  the 
increased  wavespeed  of  order  unity,  all  coincide  with  the  compressible  Rayleigh  scales  :  see 
also  section  5. 

Finally  here,  a  combination  of  the  small- r  features  in  (b)  above  and  those  for  large  r  in 
(3.14)  indicates  the  existence  of  an  0(1)  value  of  the  effective  surface  temperature  at  which 
fi-neutral  is  a  minimum,  for  the  supersonic  2D  and  less  inclined  3D  waves.  This  is  as  shown 
in  the  figures.  The  corresponding  0  value  in  the  2D  case  for  instance  represents  a  critical 
value  of  the  scaled  frequency,  below  which  2D  modes  remain  stable,  in  supersonic  flow,  and 
above  which  they  are  unstable  for  a  range  of  surface  temperatures. 

4.  THE  SEVERE-COOLING  REGIME 

This  regime  has  much  reduced  surface  temperatures  and  its  structure  is  inferred  from  the 
properties  found  in  section  3(b).  Thus  for  smail  r  the  characteristic  thickness  of  the  outer 
tier  deceases  like  r6/5,  as  does  the  streamwise  length  scale,  while  the  buffer  thickness  which 
is  proportional  to  T 2  decreases  like  r  in  view  of  (2.12a).  Along  with  this,  the  representative 
viscous-tier  thickness  falls  as  a  1^3TW  b  ecause  of  the  scaling  in  (2.7)  and  the  definition  of 
£o-  The  unsealed  thicknesses  are  therefore  of  the  orders  Re_1//2r6/5,  Re~2!3r,  Re~3^4r9^10 
respectively. 

The  severe-cooling  regime  arises  when  any  two  of  the  last  three  thicknesses  mentioned  be¬ 
come  coincident.  For  instance  if  the  buffer  thickness  coincides  with  the  viscous-tier  thickness 
then  r  ~  Re”5/6,  giving 

Tw^Re~'/2  (4.1) 
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as  the  typical  surface  temperature  here.  At  that  stage  the  other  thickness  above  also  becomes 
comparable  with  the  other  two,  and  so  a  single-tier  structure  is  expected  to  apply  for  severe 
surface  cooling.  In  effect  the  original  main  tier  and  buffer  layer  now  enter  the  viscous 
sublayer.  The  corresponding  streamwise  length  scale  is  0(Re~ 3^2)  which  is  relatively  tiny 
and  coincides  with  the  inner-tier  thickness  where  y  ~  Re~ 3'2  and  viscous  compressible  effects 
are  felt.  The  corresponding  velocity,  pressure,  temperature,  density  and  viscosity  scales  can 
be  worked  out  likewise.  The  resultant  governing  equations  for  the  severely-cooled  flow  are 
in  effect  the  unsteady  compressible  Navier-Stokes  equations, 


Pi  +  {pu)x  +  {pv)y  +  {pw)z  =  0,  (4.2a) 

r  ,  2 
p[dt  +  udx  +  vdy  +  wdz\u  =  —px  +  dx[2pux  +  (g.  -  ~/i)V.F]  +  dy[p(uy  +  ux)j  +  dz[p{wx  ■+■  uz)], 

2  (4'26) 
p[dt  +  udx+vdv+ wdz\v  =  -PyJrdy[2pvy+(p' -  -p)V .V_\+dz[p(vz  +  wy)}+dx[p.(uy+vx)]}  (4.2 c) 

2 

p[dt  +  udx  +  vdy  +  wdz]w  =  —pz  T  dz[2pw;  -f  (g!  -  -p)V.V_]  +  dx{p{wx  +  uz)]  -f  dy[p(vz  +  tav)], 

(4. 2d) 

p[Ht  +  ullx  4-  vHy  -f  wHz]  =  V.(k  grad  T  +  tu),  (4.2e) 

H  =  -(u2  ■+ v7  +  w2) p~CT,  p  =  7lpT,  (4.2f,g,h) 

/  p 

where  H  is  the  enthalpy,  r  is  the  stress  tensor  and  a  is  the  Prandtl  number.  Here  again 

the  nonlinear  interaction  is  set  up  first,  with  small  perturbations  then  being  governed  by 

the  linearized  version.  The  basic  undisturbed  flow  has  u  =  T’u,c2cJ"1[(  1  +  2 cxy)112  -  1]  and 
T  =  7^(1  +  2c13/)1/'2,  as  in  (2.3)  in  effect,  and  the  boundary  conditions  on  (4.2)  include 


U  =  V  =  w  =  0,  T  -  Tw,  at  y  -  0,  (4.3a) 

u  ~  Tw—[(  1  +  2cxy)x/2  —  1],  T  ~  7^(1  +  2 Cjy)1/2,  p  — >  0,  as  y  — >  oo,  (4.36) 

C\ 

where  Tw  is  the  scaled  surface  temperature. 

The  main  difference  from  the  previous  cooling  regime  is  that  the  viscous  Navier-Stokes 
motion  is  now  fully  compressible  (and  single-structured),  contrasting  with  the  quasi-incom- 
pressible  sublayer  properties  holding  before.  The  scales  involved  are  extreme  but  they  can 
be  verified  from  an  order-of-magnitude  argument.  Further,  the  present  flow  structure  applies 
for  any  Mach  number,  in  line  with  the  closeness  of  the  Navier-Stokes  structure  to  the  surface; 
and  both  the  spatial  and  temporal  growth  rates  induced  for  small  perturbations  are  an  order- 
of-magnitude  larger  than  those  found  previously  for  inviscid  modes,  because  of  the  tiny  scales 
involved  here. 
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This  new  predominance  of  the  viscous  modes  is  due  physically  to  the  high  surface  shear 
stresses  and  surface  heat  transfer  rates  that  the  cooling  of  the  surface  provokes.  The 
same  comment  applies  to  section  3(b).  The  precise  behavior  is  of  course  dependent  on 
the  temperature-viscosity  law  in  operation  (unlike  the  law-independent  triple-deck  case  of 
(2.2)),  but  the  same  qualitative  behavior  seems  likely  to  hold  in  general,  for  the  same  physical 
reasons.  Linearized  or  nonlinear  solutions  in  this  regime  could  therefore  be  of  some  interest, 
as  Tw  is  varied. 

5.  ON  INVISCID  MODES  ON  COOLED  SURFACES 

The  pure  inviscid  stability  of  compressible  boundary  layers  is  controlled  by  the  compress¬ 
ible  Rayleigh  equation 

pT7' 

p"  -  -^P  -  [«2(1  -  M2)  +  /32]p  =  0,  (5.1a) 

for  the  pressure  perturbation,  with  the  length  scales  as  for  the  main  tier  in  section  2  but 
now  the  wavespeed  c  is  0(1),  so  that 

=  (Uo  —  ^MooTo-1/2,  (5.16) 

cf.  (2.15a).  The  boundary  conditions  are  to  match  with  the  undisturbed  incident  free  stream 
and  to  satisfy  the  tangential-flow  surface  condition, 

p  — >  0  [or  outgoing  waves]  as  y  — »  oo,  (5.1c) 

p  =  0  at  y  =  0,  (5. Id) 

respectively.  Computational  solutions  of  (5.1)  are  given  by  Mack,  Malik  and  others  (see 
references  in  §1)  including  some  investigations  for  cooled  surfaces.  The  investigations  trace 
the  effects  of  gradual  cooling  on  the  originally  un-cooled  inviscid  modes  and  find  some  inter¬ 
esting  features.  There  is  no  significant  change  of  scale  however.  The  present  investigation 
focuses  mostly  on  new  inviscid  modes  involving  substantial  scale  changes  and.  in  particular, 
much  increased  growth  rates. 

The  main  factor  for  cooled  surfaces  where  Tw  is  small  is  the  development  of  the  0(7'2,) 
sublayer  described  by  (2.3),  associated  with  high  wall  shear  and,  temperature  gradient.  The 
previous  investigations  mentioned  above  neglect  the  sublayer  m.  a  sense  and  seek  solutions 


of  (5.1)  on  the  scale  y  of  0(1 
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,  and  by 

(5.2a) 


the  balance  of  terms  in  (5.1a).  So  M  is  large  in  the  sublayer  and  (5.1a)  reduces  to 


,  P  k  P  __  n 

P  ^  +  UJ2  -  0> 

ly  y 


where  y  —  y  +  (2ci)  1 ,  the  parameter  k  —  acM00(2c\)  l^T^/2  is  0(1)  and  ci  is  the  heat- 

transfer  factor  defined  in  (2.3b).  Here  (5.2a)  becomes  Airy’s  equation  in  terms  of  the  variable 
-1/2 

s  =  y  and  so  the  constraints  (5.1c, d)  lead  to  the  eigenrelation 


— 2iAi(  — (4fc2)1/3)  +  5i(— (4/c2)1/3)  =  0 


(5.2  b) 


for  k.  This  has  no  real  roots.  All  the  (complex)  roots  give  damped  waves  in  fact,  but  this 
type  of  mode  serves  to  suggest  the  new  scales  for  the  second  type  of  surface-cooled  mode. 
Here  c  reduces  to  0(T^2),  and  a  ~  T~2,  as  implied  by  the  first  type,  and  this  brings  in  a 
new  effect  since  now  M  is  of  order  unity.  Thus  if  c  =  T^J2c.  and  a  =  T~2a  then  M  cx  y  ^ 
to  leading  order  and  (5.1a)  becomes 


^  f  .  ,2,  L 

V  H - -  +  Ck  = : - 

2y  (2  cry)1/2 


-  1  }p  =  0. 


(5.3a) 


A  1/2 

For  a  real,  the  substitution  y  =  c2  A/^)/(8c1)1/,2  4-  s/( 2a1/2)  converts  (5.3a)  to  the  parabolic 
cylinder  equation 

Pii  -  ( J  ~0)P=  °>  i5-3b) 

where  0  —  ac4M^3/(8c1)  and  the  boundary  conditions  are  p(oo)  =  0,dp/d3  =  0  at  s  — 
(2  —  c2A/^,)d1/2/(2c1)1^2.  The  solution  is  therefore  the  parabolic  cylinder  function  p  =  £/(a,  s) 
with  a  —  —  0  (see  Abramowitz  &  Stegun  (1964)),  for  the  outer  constraint,  while  the  inner 
constraint  requires  that  s  =  xn  where  xn  are  the  roots  of 


U  (a,  xn)  =  0. 


(5.3  c) 


The  above  reference  shows  that  there  do  exist  some  roots  with  a  being  negative,  0  positive, 
as  required,  occurring  for  positive  i„.  For  those  roots  the  values  of  a,c  then  follow  from  the 
relations 

(2-c2M2OQ)al/2/(2c1)l/2  =xn,  (5.4a) 

acAMl/(2cl)  =  A0.  (5.46) 

The  resulting  wvive  speeds  c  are  all  real,  given  by 

al'2=(/3l"  +  i„/2)(2c,)1'2,  M„c  =  2(2  n./j1'2)'1'1,  (5.5a, 4) 
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where  the  first  equation  fixes  a,  the  second  c.  Hence  this  new  type,  which  is  a  temperature- 
driven  mode  because  of  the  form  of  M,  produces  new  inviscid  neutral  modes. 

The  third  case  has  c  decreased  still  further,  to  0(TW),  which  reinstates  the  influence  of 
the  velocity  profile  JJq  as  well  as  the  temperature  at  leading  order.  For,  if  c  =  Twc,  a  —  T~2a, 
then  M  ~  T^J2  M  is  small  and  the  controlling  equation  becomes 


2M  -V  2-  n 
-J—V  ~  ai 2p  =  0, 
M 


(5.6a) 


M  =  l(Vl1/2  -  1  )2 A i  -  (5.6 6) 

(Here  yl  =  2 cqy,  a  =  2c:aj,  2 =  c2).  This  can  be  recast,  by  setting  y^2  =  p(  1  +  (),  p  = 
1  +  c / 2 A ! ,  in  the  form 

(Pa  -  2P(  -  <t>2({(  +  (5.6c) 


p(oo)  =  0,  p^  =  0  at  (  =  p  1  —  1,  (5.6d) 

where  =  2aip2  acts  as  the  eigenvalue.  This  third  type  can  be  shown  to  match  with  the 
second  type,  at  large  c,  pi,  and  similar  matching  occurs  between  the  first  and  second  types. 
Hence  again  in  this  regime  neutral  modes  are  produced. 

The  fourth  type  of  cooled-surface  inviscid  mode  to  consider  has  a  structure  suggested  by 
the  work  of  section  3(c).  Here  c  is  reduced  to  0(T2)  but  a  is  0(1),  and  the  solution  of  (5.1) 
takes  on  a  three-tier  form.  In  the  bulk  of  the  boundary  layer  where  y  is  0(1),  therefore, 
(5.1a)  applies  but  with  c  replaced  by  zero  in  the  definition  of  M  in  (5.1b),  i.e.  (2.15a)  holds 

l/2 

again  with  A/0  =  UoMoo/Tq  ,  and  now 

p  ~  P  —  -a.2 A**!]3/2  as  y  — >  0+,  (5-7) 

u 

where  the  relation  between  the  constants  A**  and  P  has  to  be  found.  The  next  tier,  closer 
to  the  surface,  has  y  —  T2y  along  with  (2.3),  so  that  to  leading  order 

T7  _  rr  1/2  kV  1  .  iC.r  _  c2[(l  T  2ciy)  !  —  ljA^oo  .  . 


M  =  TlJ2M{y)  :  M 


ci(l  +  2  Cly)V4 


(5.8a) 


P  =  P  +  T3pt  +  •  •  • 

and  consequently  the  first  two  terms  in  (5.1a)  dominate,  giving  on  integration 

d]h  =  _  aM**c?/2  •  2 

dv  2  '!2M2c2  ' 


(5,86) 


(5.ftc) 


2fi 


where  the  coefficient  on  the  right  side  is  chosen  to  match  with  (5.7)  and  we  normalize  such 
that  pi(0)  =  0.  Hence  dpi/dy  decreases  like  y3  at  small  y,  which  leads  on  to  the  lowest  tier 
even  closer  to  the  surface  where  y  =  T3y  with  y  of  0(1).  There 


V  =  P  +  T*p  i  +  •  •  • , 

and  the  small  wavespced  c  =  T2c  takes  effect  since 

M  =  T3/2M  +  •  •  •  with  M  =  ( c2y  -  cjAfoo, 


(5.9a) 


(5.96) 


from  (5.1b)  with  (2.3a).  Hence  the  governing  equation,  from  (5.1a),  is 

d2pi  2  dM  dpi  2 

- - - =  CL  P. 

dy2  M  dy  dy 

The  solution  for  dpi/dy  merging  with  (5.8b,c)  at  large  y  is  then 


dpi 

dy 


a2  P  0(2J4*’r3'/2 

a  r  /  -  -s  a  n  C1  /  -  -\2 

-{c2y  -  c) - 21/2c2  [c2y  ~  c )  , 


c2 


and  so  the  required  surface  condition  (5. Id)  is  satisfied  at  y  —  0  if 

a-  -  2^p 

A  ~  3/2-  ■r' 

Ci  C 


(5.9c) 


(5.9  d) 


(5.10) 


This  provides  the  relation  for  A**,  to  be  substituted  into  (5.7).  The  resulting  problem  for  p, 
that  is,  (2.15a,b)  with  (5.7),  (5.10),  is  equivalent  to  the  one  derived  and  solved  numerically 
in  section  3(c):  see  (3.13b),  where  the  parameter  fa3  involved  for  the  two-dimensional  case 
is  in  agreement  with  that  inferred  from  (5.7),  (5.10),  with  account  taken  of  the  different 
temporal  scaling  in  (2.8).  Hence  the  stability  properties  in  subsonic  and  supersonic  flow,  for 
this  fourth  type,  are  exactly  as  given  in  Figure  4  with  the  appropriate  change  of  notation. 
Like  all  the  types  described  in  this  section,  this  new  type  appears  to  be  quite  distinct  from 
any  cooled  modes  found  previously. 


6.  FURTHER  COMMENTS 

The  main  conclusion  from  the  present  theoretical  research  [in  particular,  §§3,5]  is  that 
moderate  surface  cooling  can  enhance  the  TS  viscous-mode  activity  to  the  extent  that  the 
spatial  growth  of  the  viscous  modes  can  become  comparable  with,  or  even  exceed,  that  of 
inviscid  ones,  thus  reversing  the  normal  order.  In  addition,  the  cooling  can  completely  desta¬ 
bilize  otherwise  stable  viscous  modes  and  inviscid  modes,  at  any  Mach  number,  including  the 
production  of  radiating  waves  in  the  supersonic  realm.  In  those  two  senses,  then,  cooling  of 
the  surface  seems  very  “dangerous”  with  regard  to  transition  in  the  compressible  boundary 
layer. 
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The  TS  growth,  when  it  exceeds  that  of  the  pure  inviscid  modes  studied  previously, 
occurs  for  disturbances  which  are  confined  relatively  close  to  the  surface,  as  opposed  to  the 
full  inviscid  ones  stretching  all  the  way  across  the  boundary  layer.  In  particular,  for  all 
the  new  modes,  the  high-heat-transfer  sublayer  plays  a  substantial  active  part  in  enhancing 
the  viscous  growth,  through  the  correspondingly  high  shear-stress  production  at  the  surface. 
It  could  be  for  that  reason,  i.e.  the  appearance  of  multi-structured  responses  within  the 
boundary  layer,  that  previous  theoretical,  computational  and  experimental  studies  of  surface¬ 
cooling  effects  appear  to  have  not  detected  the  new  mode  properties,  although  there  are  other 
explanations  which  suggest  themselves  such  as  considerable  finite-Reynolds-number  effects; 
see  also  below.  Experiments  and  computations  focusing  on  the  possibility  of  such  strong  near¬ 
surface  behavior  remain  to  be  done,  it  seems,  and  they  could  be  of  much  interest  scientifically 
and  for  practical  applications.  There  is  some  measure  of  agreement,  nevertheless,  between 
the  present  theory  and  the  Lysenko  and  Maslov  (1984)  experiments.  Thus  the  theory  predicts 
that  as  the  temperature  ratio  Tw  falls  from  the  non-cooled  regime  to  the  moderate-cooling 
regime,  for  first  modes  at  fixed  Reynolds  number,  the  typical  streamwise  wavenumber,  aD 
say,  and  the  typical  frequency,  uD  say,  increase  as  T~ 3^2,  Tj1  respectively:  see  §2.  So  for  the 
1  to  0.76  decrease  in  Tw  examined  in  Lysenko  and  Maslov’s  Figure  2  the  predicted  increase  in 
uj°  is  by  a  factor  1.316,  which  is  not  inconsistent  with  the  measured  increase  in  the  neutral 
value  of  U)D  shown  in  that  figure.  (See  also  comparisons  for  the  insulated  case  in  Smith 
1987).  Similar  trends  for  decreasing  Tw  are  observed  in  their  Figure  5.  Again,  Lysenko  and 
Maslov’s  Figure  4  has  a  case  of  a  1  to  0.79  decrease  in  Tw,  for  which  the  present  theory 
implies  an  increase  by  a  factor  (0.79)"3/2,  i.e.  1.398,  in  the  corresponding  spatial  growth 
rates  —  This  also  is  not  inconsistent  with  the  measured  growth  curves  in  their  Figure  4. 
Hence  as  far  as  their  Figure  4  is  concerned  there  appears  to  be  some  agreement  both  on  the 
delay  of  instability  and  on  its  increased  growth  rates,  with  cooling  of  the  surface,  although 
their  interpretation  (see  their  figures)  is  that  the  growth  rates  in  fact  decrease.  The  latter 
would  seem  to  be  due  to  extra  normalization  factors  used,  whereas  the  theory  indicates  that 
the  effects  of  the  temperature  ratio  Tw  should  be  the  overwhelming  ones  eventually,  as  Tw 
falls. 

A  further  note  of  caution  should  be  added  here,  concerning  nonparallel-flow  effects.  We 
recall  that,  strictly,  supersonic  first  modes  within  the  wave-Mach  cone  are  stable  and  gov¬ 
erned  at  leading  order  by  nonparallel-flow  effects,  for  Tw  of  order  unity  (Smith  1987).  There¬ 
fore  the  experimental  results  (of  Lysenko  and  Maslov  (1984),  and  others)  inside  the  wave- 
Mach  cone  could  well  be  influenced  to  a  significant  extent  by  nonparallelism  rather  than  being 
true  instability  measurements,  and  they  may  in  turn  miss  the  new  supersonic  instabilities 
predicted  for  cooled  surfaces. 
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The  corresponding  nonlinear  processes  at  higher  disturbance  amplitudes  which  are  noted 
earlier  in  the  text  should  also  be  of  much  interest  with  respect  to  transition  in  the  cooled 
boundary  layer.  For  example,  the  nonlinear  break-ups  described  by  Smith  (1988)  then  be¬ 
come  attainable  within  finite  time;  and  forms  of  vortex-wave  interaction  can  also  be  induced, 
including  those  of  Smith  &  Walton  (1989)  and  Hall  &  Smith  (1989).  As  well  as  the  nonlinear 
problems  mentioned  in  the  text,  that  inferred  from  section  3(c),  for  reduced  surface  tem¬ 
peratures,  is  governed  by  the  three-dimensional  unsteady  boundary-layer  equations  coupled 
with  the  pressure-displacement  law  which  is  given  in  transformed  terms  by  (3.5b).  This  can 
be  expressed  as  an  integral  relation  between  the  unknown  pressure  and  displacement,  in  two- 
dimensions  or  three-dimensions.  Further,  the  finite  time  for  nonlinear  break-up  referred  to 
above  becomes  shorter  as  the  surface  is  cooled,  suggesting  that  intermittency  (see  references 
above)  may  arise  that  much  sooner  on  a  cooled  surface. 

There  are  many  other  regimes  of  concern.  These  include  the  hypersonic  realm  with  its 
strong  vorticity  mode  (Brown  &  Smith  1989);  the  transonic  realm  and  its  peculiar  instability 
features  (Bowles  &  Smith  1989);  cross-flow  effects;  upper-branch  properties,  e.g.  between 
the  mounds  of  instability  discussed  in  section  3(c);  other  temperature-viscosity  laws,  given 
that  the  Chapman  assumption  is  vital  in  the  new  near-surface  layers,  although  again  other 
laws  of  concern  seem  likely  to  induce  broadly  similar  behavior;  and  the  upstream-influence 
analogue  which  is  currently  being  considered.  These  additional  aspects  seem  less  significant 
in  the  long  run,  however,  than  the  issues  associated  with  nonlinearity,  given  the  increased- 
instability  properties  found  here  for  cooled-surface  boundary  layers. 
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2/3  buffer  tier  and  a  3/4  viscous  lower  tier.  Further  cooling  shrinks  the  instability  structure 
still  further,  as  described  in  sections  3(b)  and  4. 


Figure  2(b). 


Figure  2(b).  Neutral  value  of  r  versus  frequency  fi  for  M »  =  2  and  f3/a  —  0,1, 1.25, 1.5, 1.75. 
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Figure  2(c).  Neutral  value  of  a  versus  frequercy  fi  for  M 


Figure  2(f)- 
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Figure  2(f).  Growth  rate  versus  frequency  fl  for  M0 0  —  0.5,  /3  =  0  and  r  —  0.4, 0.5,  1. 
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.  Asymp 
and  t  = 


Figure  3(d) 
for  8  =  0  < 


Figure  4(a),  Asymptotic  solution  for  Re(a )  as  a  function  of  f  =  r/fi,  for  r  >>  1  and 
>>  1  for  (3  —  0  and  M*  =  2,  for  r  ~  0(1),  Q  >>  1  and  (3  =  0  from  (3.12),  along  with 
the  full  computational  solution  for  M <*,  =  2,  (3  =  0  and  r  =  10. 
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Figure  4(f)- 


Figure  4(f).  Neutral  value  of  f  versus  frequency  Q,  for  M =  2  and  (3 /a  =  1. 
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Figure  4(g). 
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Figure  4(g).  Neutral  value  of  f  versus  frequency  Q  for  =  2  and  /3/a  =  1.5. 
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16  Abstract  inf^uence  0f  surface  cooling  on  compressible  boundary-layer  instability 

is  discussed  theoretically  for  both  viscous  and  inviscid  modes,  at  high  Reynolds 
numbers.  The  cooling  enhances  the  surface  heat  transfer  and  shear  stress,  creating 
a  high-heat-transfer  sublayer.  This  has  the  effect  of  distorting  and  accentuating 
the  viscous  Tollmien-Schlichting  modes  to  such  an  extent  that  their  spatial  growth 
rates  become  comparable  with,  and  can  even  exceed,  the  growth  rates  of  inviscid 
modes,  including  those  found  previously.  This  is  for  moderate  cooling,  and  it  ap¬ 
plies  at  any  Mach  number.  In  addition,  the  moderate  cooling  destabilizes  otherwise 
stable  viscous  or  inviscid  modes,  in  particular  triggering  outward-traveling  waves 
at  the  edge  of  the  boundary  layer  in  the  supersonic  regime.  Severe  cooling  is  also 
discussed  as  it  brings  compressible  dynamics  directly  into  play  within  the  viscous 
sublayer.  All  the  new  cooled  modes  found  involve  the  heat-transfer  sublayer  quite 
actively,  and  they  are  often  multi-structured  in  form  and  may  be  distinct  from  those 
observed  in  previous  computational  and  experimental  investigations.  The  correspond¬ 
ing  nonlinear  processes  are  also  pointed  out  with  regard  to  transition  in  the  cooled 
compressible  boundary  layer.  Finally,  comparisons  with  Lysenko  and  Maslov's  (1984) 
experiments  on  surface  cooling  are  presented.  / 
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